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Abstract. We found that the order for the crumpling transition of an intrinsic 
curvature model changes depending on the distance between two boundary vertices 
fixed on the surface of spherical topology. The model is a curvature one governed 
by an intrinsic curvature energy, which is defined on triangulated surfaces. It was 
already reported that the model undergoes a first-order crumpling transition without 
the boundary conditions on the surface. However, the dependence of the transition 
on such boundary condition is yet to be studied. We have studied in this paper this 
problem by using the Monte Carlo simulations on surfaces up to a size N = 8412. The 
first-order transition changes to a second-order one if the distance increases. 
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1. Introduction 

Over the past two decades, a considerable number of studies have been performed on 
the phase structure of the surface model of Helfrich 1J, Polyakovj2], and Kleinert j3] 

for biological membranes HSEEISISIiniimilSIISIll!. The surface models can 
be classified into two groups, which are characterized by the curvature energy in the 
Hamiltonian; one is an extrinsic curvature model [TBllTo1IT?llT^ITDll2T)l I2T] 1221 1231 12311251 
E3 E3 EB1 EH1 E01 EJJ E2], and the other is an intrinsic curvature model [SSI El E51 EE] ■ 

The extrinsic curvature model is known to undergo a first-order transition between the 
smooth phase and the crumpled phase on tethered spherical surfaces [22J ESI El] • 

Studies have also focused on the phase structure of the model with intrinsic 
curvature [37J EH ESI HOI EJJ- It has been shown that the first-order transition can 
be seen in spherical fluid/tethered surfaces [SHIES], tethered surface of disk topology 
[ID] , and tethered surface with torus topology [UJ. 

However, little attention has been given to the boundary conditions of the surface 
models, although the center of the surface can be fixed to prevent translation in Monte 
Carlo (MC) simulations. Though, it is as yet unclear how boundary conditions of the 
surface influence the crumpling transition. 

Therefore, we should carefully study the influences of boundary conditions on the 
phase transition. In fact, we know that a phase transition is significantly influenced by 
boundary conditions, which fix some of the dynamical variables to a prescribed value. 
Moreover, it is clear that artificial vesicles can be supported with substrates such as 
glass plates or beads in an aqueous solution. 

In this paper, we study how a boundary condition influences the crumpling 
transition of the tethered surface model with intrinsic curvature on spherical surfaces. 
The boundary condition is imposed on the surface with two fixed vertices of distance 
2L(N), which depends on the total number of vertices iV of the surface. We will find that 
the order of the transition changes from first-order to second-order and higher-orders 
in the limit of N — > oo when L(N) is increased from Lq(N) to 1.5Lq(N) and 2Lq(N), 
respectively, where L (N) is a radius of the initial sphere such that the Gaussian energy 
{Si)/N is approximately equal to (Si)/N ~ 3/2 at the start of MC simulations. The 
result in this paper is in sharp contrast to that of the fluid surface model with extrinsic 
curvature, where the phase transition is strengthened with the increasing L(N) [HJJE2]- 

2. Model and Monte Carlo technique 

The tethered surface model is defined by the partition function 



where a is the curvature coefficient, and /' denotes the boundary condition in which 
two vertices are fixed and separated by a distance of 2L(N). S(X) denotes that the 




»=i 



(1) 



S{X) = S 1 + aS 3 
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Hamiltonian S depends on the position variables X of vertices. Si and S3 are defined 
by 

S 1 = J2(X l -X J )\ S 3 = -EMV27r), (2) 

{») » 

where Yluj) i n Si is the sum over bonds (ij) connecting the vertices i and j. The bonds 
(ij) are the edges of the triangles. 5i in Eq. (j2j) is the vertex angle, which is the sum 
of the angles meeting at the vertex i. We call S 3 the deficit angle term, because <5j — 2n 
is just the deficit angle. We note that 27r) is constant on surfaces of fixed genus 

because of the Gauss-Bonnet theorem. 

We comment on the unit of physical quantities in the model. By letting a be a 
length unit in the model, we can express all quantities with unit of length in terms of 
a. Hence, the unit of Si is a 2 , and that of L Q (N) is a. We fix the value of a to a — 1 in 
this paper, because the length unit can be arbitrarily chosen in the model. The unit of 
a is expressed by kT, where k is the Boltzmann constant, T is the temperature. Note 
also that varying the temperature T is effectively identical to varying a in the model. 

Triangulated surfaces are obtained by dividing the icosahedron. By splitting the 
edges of the icosahedron into ^-pieces, we have a surface of size N = 10£ 2 + 2. These 
surfaces are identical to those in |2Sl 121! • These surfaces are characterized by N 5 = 12 
and iV 6 = iV— 12, where N q is the total number of vertices with a co-ordination number 

The radius of the sphere is fixed to Lq(N), which depends on N and was defined 
so that the mean Gaussian energy (Si)/N is approximately equal to (Si)/iV~3/2 as 
mentioned in the Introduction. We assumed the following values of Lq(N): Lq(N) = 14 
for N = 1442, L (N) = 18.6 for N = 2562, L (N) = 25.6 for N = 4842, and L (N) = 34 
for iV = 8412. Under those values of Lq(N) for the radius, the relation Si/N = 1.5 is 
almost satisfied in the initial configurations for MC simulations. Note that the value 
of the radius L (N), which satisfies (Si)/N ~ 3/2, depends on the construction of the 
surface. Note also that the sphere with the radius Lq(N) seems correspond to a real 
physical membrane at sufficiently low temperature or in the limit of a — > 00. 

The distance 2L(N) between two vertices is fixed to three different values such that 

L(N) = L (N), 1.5Lo(N), 2L (N) (3) 

where L (N) is a radius of an initial sphere constructed from the icosahedron as 
described above. Then, the distance between two fixed vertices on the surface with 
L(N) = 2L (N) is just 4Z/ (iV) for example. The boundary conditions are thus imposed 
on the model by these two fixed vertices of distance 2L(N). 

If it were not for the boundary condition defined in Eq. (jH)), we have (Si)/N — 3/2, 
which comes from the scale invariant property of Z. However, as we will see later, this 
relation is slightly broken due to the boundary condition. 

Figure ^a) shows a schematic view of the diameter 2L Q (N) of a sphere, and Fig. 
^b) shows that of the distance 2L(N) between the two vertices of the expanded spherical 
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2L (N) 2L(N) 
(a) (b) 

Figure 1. Schematic drawing of (a) the diameter 2Lq(N) of a sphere which is the 
starting configuration for the MC simulations, and that of (b) the distance 2L(N) 
between two vertices of the expanded spherical surface, where L(N) is fixed to three 
different values L(N) = L (N), L(N) = l.5L {N), and L(N) = 2L (N) in the 
simulations. Surfaces were drawn as simple as possible, since they are aimed only 
at showing Lq(N) and L(N). The symbols (•) denote two fixed vertices on each 
surface. The surface in (b) is obtained by expanding the distance of the vertices from 
2L (N) to 2L(N) during thermalization in MC. 



surface. Surfaces were drawn as simple as possible, since they are aimed only at showing 
L Q (N) and L(N). 

The procedure for choosing and fixing the position of the boundary vertices is as 
follows: Firstly, we choose a pair of vertices which are on a straight line passing through 
the center of the triangulated sphere for the starting configuration. The canonical x- 
coordinate axis in R 3 = {(x, y, z)\x, y, z G R} is chosen as the straight line in the 
surface. The distance between these two vertices on the axis is 2Lq(N) in the starting 
configuration. Secondly, the distance is expanded from 2L (N) to 2L(N) along the axis 
during the thermalization in MC, where L(N) is given in Eq. (JHJ). The total number of 
the thermalization in MC is about 1 x 10 7 for expanding the distance. Some additional 
thermalization MCS are performed after the expansion, if necessary. 

When the condition L(N) = 2Lq(N) was chosen, the distance 2L(N) is expanded 
from 2L (N) = 14 to 2L(N) = 28 during the thermalization for the N = 1442 surface. 
The two fixed vertices of distance 2Lq{N) = 14 are shifted by distance 7 x 10~ 5 at 
every 100 MCS to the opposite direction with each other along the axis during the first 
1 x 10 7 thermalization MCS. Note also that the results of the simulation is completely 
independent of how the surface with the distance 2L(N) is constructed from the starting 
configuration with the distance 2L (N) between the boundary points. 

We comment on a relation between the deficit angle term S3 and the integration 
measure dXi of the partition function in Eq. (JI}. The integration measure fli dXi can 
be replaced by FJi Q.fdX^ where is the co-ordination number of the vertex i jl2]. This 
a is believed to be 2a = 3. Then, it is possible to consider that qf is the volume weight 
of the vertex i in the measure cOQ. Hence, we can extend 2a to non-integer numbers 
by assuming that the weight is arbitrarily chosen. Moreover, we can also extend qi to 
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continuous numbers for the same reason. Hence, the weight FL it can De replaced by 
Til &i =exp(a Yji log^i). Including a constant weight (27r) a , we have S 3 in Eq. (J2J). 

It should also be noted that S3 (5) in Eq. (j2J can make the surface smooth not only 
in the model with the Gaussian term but also in the Nambu-Goto surface model within 
the class of tethered surfaces [HE]- Whereas the term 5*3 (g) = — Yh l°g<?i is constant on 
the tethered surfaces, and moreover the equivalent term S 3 (g) = J2i(li~^) 2 plays no role 
in smoothing fluid surfaces [37]. We therefore have discovered a significant difference 
for the role of S3 (5) in Eq. (J2J) and that of the corresponding term Ss(q) for smoothing 
the surface. 

The vertices X are shifted so that X' = X+5X, where 5X is chosen randomly in a 
small sphere. The new position X' is accepted with the probability Min[l, exp(— AS)], 
where AS = S(new)— S(old). We use two sequences of random numbers called Mersenne 
Twister |33] ; one for three dimensional random shift of X and the other for Metropolis 
accept /reject. The radius of the small sphere for the shift SX is chosen so that the rate 
of acceptance for X is about 50%. We introduce the lower bound 1 x 10 7 for the area 
of triangles. No lower bound is imposed on the bond length. 



3. Results 

Monte Carlo simulations were performed on the surfaces of size N = 1442, iV = 2562, 
N = 4842, and N = 8412. The total number of sweeps of MC simulations were about 
0.8 x 10 8 ~ 1.2 x 10 8 at transition point a c for each surface. Relatively small numbers 
of sweeps were performed at non-transition points a=fot c . 
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Figure 2. X 2 vs. a obtained under the conditions (a) L(N) = Lq(N), (b) 
L(N) = 1.5L (N), and (c) L(N) = 2L (N). Solid lines were obtained by multihistogram 
reweighting technique. 



First, we plot the mean square size X 2 in Figs. Efa)H2tc), where X 2 is defined by 
X 2 = ^Y.{X t -X)\ X = Wx„ (4) 

i i 

where X is the center of the surface. X 2 reflects the size of surfaces even on surfaces with 
the conditions given by Eq. (J3J) for two fixed vertices. The boundary conditions for these 
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data in Figs. H», Mb), andEfc) are given by L(N) = L (N), 1.5L (N), and 2L (N), 
respectively. Solid lines drawn on the data were obtained by multihistogram reweighting 
technique |44j . The reweighting analysis was done by using all the simulation data in 
Fig. Ufa), and it was done by using some of the data in Figs. Efb)|Hc). 

FigureEfa) shows discontinuous changes of X 2 against a for L(N) — L Q (N), whereas 
X 2 continuously changes in Figs. I2fb)|2tc) corresponding to L(N) = 1.5L (iV) and 
L(N) =2L (N). The discontinuity of X 2 indicates that the model undergoes a first- 
order transition under the condition L(N) = Lq(N). 
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Figure 3. Si/N vs. a obtained under the conditions (a) L(N) = Lq{N), (b) 
L(N) = l.5L (N), and (c) L(N) = 2L (N). Solid lines were obtained by multihistogram 
reweighting technique. 



We plot the Gaussian energy Si/N against a in Figs. Et a )HSf c )- As mentioned 
in the previous section, S\/N should be equal to 3/2 whenever no specific boundary 
conditions are imposed on the model. However, Si/N may deviate from 3/2 in Fig. El 
because the surface is fixed with two vertices separated by a distance 2L(N). In fact, 
we see a discontinuous change of Si/N in Fig. Efa), which corresponds to the condition 
L(N) = L (N). The discontinuity of Si/N indicates that the model undergoes a first- 
order transition under L(N) =Lq(N). We can also find that Si/N deviates from 3/2 in 
Figs, lib) He). 

The bending energy S^/Nb is plotted in Figs. IHa)-|lfc), where iV^ is the total 
number of bonds; N B = 3N — 6. S 2 /N B reflects the smoothness of the surface, although 
it is not included in the Hamiltonian. A phase transition can be called a first-order one 
if some physical quantity discontinuously changes. Therefore, the discontinuous change 
of S 2 /N B in Fig. Ufa) also supports the first-order transition of the model under the 
condition L(N) = L (N). 

The specific heat Cs 2 for the bending energy S 2 is defined by 

Cs 2 = ^<(S 2 -<S 2 >) 2 >, (5) 

which is the variance of S 2 - Note that the coefficient b 2 , which is the squared bending 
rigidity, is excluded from the right-hand side of Eq. (jSJ) . The reason for this is because 
the term bS 2 is not included in the Hamiltonian. However, we expect that an anomalous 
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Figure 4. The bending energy S2/NB vs. a obtained under the conditions (a) 
L{N) = L (N), (b) L(N) = 1.5L (N), and (c) L(N) = 2L (N). S 2 is not included 
in the Hamiltonian. Nb is the total number of bonds. 
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Figure 5. The specific heat Cs 2 vs. a obtained under the conditions (a) L(N) 
L Q (N), (b) L{N) = 1.5L (N), and (c) L(N) = 2L {N). 



behavior of Cs 2 , which is typical of the phase transition, is independent of whether the 
coefficient b 2 is included in Cs 2 or not. 

Figures EJa), E(b) and EJc) show C*5 2 against a obtained under the conditions 
L(N) =Lq(N), 1.5Lq(N), and 2Lq(N), respectively. The error bars on the symbols in 
the figures are the statistical errors, which were obtained by the so-called the binning 
analysis. The solid lines were drawn by the multihistogram reweighting technique. We 
see the expected anomalous behavior of Cs 2 in Fig. Efa), which reflects a discontinuous 
transition. It is also clear from Figs. EJb) and^c) that the transition is softened with 
the increasing distance L(N); the peak values of Cs 2 become lower and lower when 
L(N) increases. 

We obtained the peak value C!? ax of the specific heat by the multihistogram 
reweighting technique. The statistical errors were also obtained by that technique and 
shown in Figs. Eta)HEtc) with the error bars. The errors were almost invisible in Fig. 
EJa), because they were relatively smaller than those shown in Figs. E^b) and|^c). 
The reason why the errors in Figs. 10(b) and E(c) are relatively large seems that the 
number of data point in the simulations is not so large; overlapping of energy is slightly 
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Figure 6. The peak values C5? ax vs. N obtained under the conditions (a) 
L(N) = Lq(N) in a log-log scale, (b) L(N) = 1.5Lq(N) in a log-log scale, and (c) 
L(N) = 2Lq(N) in a linear-log scale. C<? ax and the statistical errors were obtained by 
the multihistogram reweighting. The error bars denote the statistical errors. 



insufficient for the reweighting analysis in those cases. 
Cj 2 ffl is expected to scale according to 

C£ ax ~ N a , (6) 

where a is a critical exponent of the transition. By fitting the data in Figs. HJa) and 
Efb) to Eq. (JHJ), we have 

a! = 1.15 ± 0.10, [L(N)=L (N)), 

a 2 = 0.42 ± 0.06, [L(N) = 1.5L (N)}, (7) 
a 3 = 0.35 ± 0.07, [L{N)=2L (N)\. 

The value of o~\ = 1.15(10) indicates that the transition is of the first order, whereas 
o"2 = 0.42(6) and 03 = 0.35(7) indicate second-order transitions. 

The values <t 2 = 0.42(6) and <r 3 = 0.35(7) in Eq. (JJJ) can be compared to a known 
value a = 0.58(10) of the model with extrinsic curvature reported in ^7]. However, we 
have no definite conclusion about whether the two models are in the same universality 
class, because o~ 2 and o~ 3 slightly deviate from a = 0.58(10). Nevertheless, it is possible 
that two models are in the same class, because a of our model changes depending on 
the distance L(N). It is also expected that the transition disappear when L(N) further 
increases. We should note that a possibility of a discontinuous change of a against L(N) 
in the limit of N — > 00 is not eliminated. 

The first-order transition should be reflected in S3/N such that S3/N 
disco ntinuously changes at the transition point. In order to see this, we plot the 
distribution (or histogram) of the intrinsic curvature energy S3/N in Figs. 0(a)-Efc). 
We find that h(S 3 ) has a double peak structure in Fig. Eta.) under the condition 
L(N) = L (N). This clearly shows that S3/N is discontinuous and that the model 
undergoes a first-order transition at a = 750 on the iV = 8412 surface. On the contrary, 
S 3 /N smoothly changes in Figs. E(b)E( c ) under the conditions L(N) = 1.5L (N) and 
L(N) = 2L (N). 




Figure 7. The histogram h(S3) of S3/N obtained at the transition point of the 
N = 8412 surface under the conditions (a) L(N)=L (N), (b) L(N) = 1.5L (N), and 
(c) L(N) — 2Lq(N). A double peak structure of h(S3) in (a) indicates a first-order 
transition. 



We can also see a discontinuous change of S 3 /N at a = 750 in a plot of S 3 /N 
against a, which is not depicted here. However, the discontinuity of S3/N is not so 
clear. Therefore, we can hardly clarify the order of the transition of the model under 
L(N) =Lq(N) without the histogram h(Ss) shown in Fig. Ufa). We understand that the 
transition is softened by the conditions of Eq. (JHJ) although it remains in the first-order 
transition only under the condition L(N) = L (N). 



(a) L = L Q , a = 750 (smooth) (b) L = L , a = 750 (crumpled) 




(c) The surface section of (a) (d) The surface section of (b) 

Figure 8. Snapshots of surfaces in (a) the smooth phase and (b) the crumpled phase 
at the transition point a — 750 of the TV = 8412 surface under L(N)=Lq(N), (c) surface 
sections in (a), and (d) surface sections in (b). 



Snapshots of surfaces in the smooth phase and the crumpled phase are respectively 
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shown in Figs. Ufa) andlHl^b), which are obtained at the transition point a = 750 of the 
N = 8412 surface under L(N) = L (N). The surface sections are shown in Figs. lEfc), 
Btd). 



(a) L = 1.5L , a = 480 (b) L = 2L , a = 460 





(c) The surface section of (a) (d) The surface section of (b) 

Figure 9. Snapshots of surfaces (a) at the transition point a = 480 of the N = 8412 
surface under L(N) = 1.5L (N), (b) at the transition point a = 460 of the N = 8412 
surface under L(N) =2Lq(N), (c) surface sections in (a), and (d) surface sections in 
(b). 

Figures E^a) and Efb) are snapshots of surfaces obtained at the transition point 
a = 480 of the iV = 8412 surface under L(N) = 1.5L (N) and at the transition point 
a = 460 of the iV = 8412 surface under L(N) = 2Lo(N), respectively. The surface sections 
in Figs. E(a) andEJb) are shown in Figs. H2c) andEfd), respectively. 

We understand from the snapshots in Figs. IHlandlHlthat the surfaces become oblong 
if the distance L(N) is increased. Therefore, it is possible to extract the string tension 
of surfaces from MC data if L(N) is sufficiently large, although the model is defined on a 
fixed connectivity surface. However, the string tension of the model will be continuous at 
the transition point because the transition becomes a second-order or higher-order one 
when L(N) is sufficiently large. This implies that the string tension does not play the 
role of an order parameter in the crumpling transition of a fixed connectivity surface 
model. However, it is very interesting that the strength of the crumpling transition 
changes depending on the boundary condition. 
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4. Summary and conclusions 

We have investigated the phase structure of a surface model with an intrinsic curvature 
under a boundary condition such that two fixed vertices are separated by a distance 
2L(N). The model was known as the one that undergoes a first-order crumpling 
transition without the boundary condition [SHI HIH ISj- This paper aimed to show 
how boundary conditions influence the phase transition, and we performed extensive 
MC simulations on the spherical tethered surfaces up to a size iV = 8412. 

The distance L(N) was chosen to be three different values for each N; L(N) = 
L {N), L(N) = 1.5Lo(JV), and L(N) = 2L (N), where L (N) depends on N and is the 
radius of the initial sphere at the start of the MC simulations. We used the following 
values of L (N): L (N) = 14 for N = 1442, L (N) = 18.6 for N = 2562, L (N) = 25.6 
for iV = 4842, and L (N) =34 for iV = 8412, so that the relation Si/N=1.5 is almost 
satisfied at the starts. 

We have found that the transition is softened as the distance L(N) is increased. 
The order of the transition remains in the first-order under L(N) = L (N), and it 
changes to the second-order when the distance is increased to L(N) = 1.5Lq(N) and 
L(N) = 2Lo(N), where the peak C<? ax of the specific heat for the bending energy scales 
according to C^ ax ~ N a ,o~ < 1 at the transition point. The first-order transition at 
L(N) = Lq(N) was confirmed with a double peak structure in the histogram h(Ss) of 
the intrinsic curvature energy S3/N, which is included in the Hamiltonian. 

The result is in sharp contrast to that of a fluid surface model with extrinsic 
curvature, where the crumpling transition is strengthened when the distance between 
two fixed vertices is increased under a specific condition [^H E2] at a sufficiently large 
L(N). In fact, the phase transition of the tethered surface model in this paper is softened 
at sufficiently large L(N) as stated above. Moreover, the phase transition is expected 
to disappear if L(N) was further increased, where the surface becomes sufficiently 
oblong. Therefore, we consider that string tension does not play the role of an order 
parameter in the crumpling transition if the model is a fixed connectivity one at least. 
Nevertheless, the fact that the strength of the crumpling transition changes depending 
on the boundary condition is very interesting, because the boundary condition seems 
accessible in real physical membranes. If the transition was observed in a biological 
membrane that is governed by the intrinsic curvature, the strength of the transition can 
be handled by fixing and expanding the surface. 

It is interesting to study the phase structure of the model on dynamically 
triangulated fluid surfaces under the same boundary condition, where the surface is 
expected to be oblong and linear if the distance L(N) is increased to a sufficiently 
large size. For future experimental studies on the crumpling transition in biological 
membranes, further numerical studies on the surface models will provide more detailed 
and helpful information about the influence of boundary conditions on the transition. 
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